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Abstract 

We investigate theoretically the coherent longitudinal and transversal spin relaxation of pho- 
toexcited electrons in quantum wells in quantized magnetic fields. We find the relaxation time for 
typical quantum well parameters between 10^ and fO^ ps. For a realistic random potential the 
relaxation process depends on the electron energy and 5- factor, demonstrating oscillations in the 
spin polarization accompanying the spin relaxation. The dependence of spin relaxation on applied 
field, and thus on the corresponding "magnetic" length, can be used to characterize the spatial 
scale of disorder in quantum wells. 
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I. INTRODUCTION 



A magnetic field applied to a two-dimensional (2D) electron system changes both the 
orbital and spin dynamics of the carriers. A crucial aspect of this dynamics is spin re- 
laxation, which arises due to spin-orbit (SO) coupling. This coupling is well-described in 
zincblende (OOl)-grown structures by a Hamiltonian Hso that is the sum of two terms, 
a Rashba Hamiltonian^'^ i^R = ckr {axPy — CyPx) / ^ and a Dresselhaus Hamiltonian^'^ 
Hb — ckd {cxPx ~ (^yPy) where q;r and are coupling constants, a are the Pauli ma- 
trices, and p|| = ipxiPy) = — — (e/c)A||. Here e is the electron charge, and Ay is a 
vector-potential of the external field; and Hd arise, respectively, due to the artificial 
macroscopic asymmetry of the structure, and due to the microscopic inversion asymmetry 
of the zincblende unit cell. The coupling constants ctR and typically range from 10~^° to 
10-9 eVcm.5 

The spin relaxation of conducting electrons is usually described using a Dyakonov-Perel'^ 
approach, where it is assumed that the orientation of the spin precession axis changes ran- 
domly through scattering by impurities. In the absence of an external magnetic field the 
spin relaxation rate is 7dp ~ {a/U) p'^r, with r the momentum relaxation time, and where 
a depends on «r and aj)- Since spin relaxation arises due to the random spatial motion of 
the electron, the orbital effect of a magnetic field, which can restrict the region over which 
an electron can sample the effect of impurities, influences the spin relaxation. At cUcT ^ 1, 
the relaxation rate decreases by a factor^ of cu'^t'^, where cUc — \e\B/mc is the cyclotron 
frequency for a magnetic field B, with m being the electron effective mass. In such a strong 
field the electron path becomes close to a circle, with scattering effects negligible. On this 
orbit the mean spin precession angle vanishes due to the fact that the electron velocity 
V||(t-|-7r/a;c) = — V[|(t), and, therefore, the randomness in the precession is suppressed. For a 
random SO coupling a non-quantizing magnetic field can, however, speed up spin relaxation 
and make the relaxation process Gaussian rather than exponential in its time dependence.^ 
The S— dependence of the relaxation rate in nonquantizing fields can demonstrate magne- 
toquantum oscillations, as shown by self-consistent Born approximation for 2D electron gas 
in a short-range potential of impurities.^ Analysis of spin relaxation in weak magnetic fields 
allows the extraction of SO coupling parameters from experimental data.^° 

Studies of spin dynamics of itinerant electrons typically assume that their motion can be 
described semiclassically. However, spin dynamics of carriers with quantized lateral motion 
is of interest both for understanding of the fundamental physics of spin transport^^ and for 
applications of nanosize systems in spintronics. An interesting example of a system with 
quantum lateral motion is 2D electron gas in a strong magnetic field, where Landau states 
must be used to represent the electrons and only a few (or a fraction) of the Landau levels is 
occupied. An analogous regime has been widely discussed for quantum dots, where electrons 
are confined by an external potential. There the SO coupling mixes states with opposite 
spins, making spin-fiip transitions accompanied by phonon emission possible. Here we 
are interested in itinerant two-dimensional electrons, where early experimental data has 
shown the suppression of spin relaxation in strong magnetic fields. The treatment of 2D 
electrons in quantizing fields is rather subtle, ^"^ with different theoretical techniques giving 
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different results even for static properties such as the density of states. The analytical 
approaches require approximations that might be not widely applicable. The situation 
becomes even more complicated for the response functions, such as those describing charge 
and spin currents, as well as relaxation processes. An analytical study by Bastard^^ of spin 
relaxation in the self-consistent Born approximation, for a short-range random potential and 
a small electron g^-f actor, where the Zeeman splitting is small in comparison with the level 
broadening, showed that both spin-orbit coupling and disorder play a role, and allowed an 
estimation of the relaxation rate. 

Here we perform a numerical study of the problem using the exact diagonalization tech- 
nique for a large but finite-size system, and show the diverse physical mechanisms that 
contribute to the relaxation process. This approach has proven its applicability in calcula- 
tions of the spin-Hall conductivity of a disordered 2D electron gas,^^ where the spectrum 
and the full set of eigenstates are required for the calculation. 



II. MODEL OF DISORDER AND SPIN DYNAMICS 



We consider a quantum well in a magnetic field parallel to the growth direction, B = 
Bz. The field forms spin-split Landau levels, inhomogeneously broadened by the disorder. 
In undoped quantum wells the random potential arises due to the monolayer islands of 
the parent compounds at the interfaces (thickness fluctuations) and due to the content 
variations near the interfaces. The island patterns depend on the growth conditions. 
The Hamiltonian for an electron in an undoped quantum well with disorder is 

H^^ + liB9{(r-B)/2 + Hso {cr, py) + U {p) , (1) 

where the random contribution U (p) = ^ (P ~ R-rf) 5 where is the position of a 
defect. As a model, we consider the Gaussian potential V (p) = exp {—p'^/Rf} with the 
areal density of defects (/?) = Xld^ (p ~ ^d) ? and a correlation function (A^ (0) A^ (p)) = 
Aimp^ (p) ! where Ai^p is the average areal density of defects. To avoid the uniform shift of 
the Landau levels, we assume that Vr = ^Vg varies from site to site, being either positive 
or negative such that the mean value {U (p)) = 0. The correlation function of the random 
potential is {U (0) U (p)) = (U^) F^{p), where^^ 

{U')-In,^,RX^ F^ip) ^ e-^y^^l (2) 

At B — the momentum relaxation time in this model is given by: 

r 27r^Rl, A>i?5, 
- - N,^.KR.^s (3) 

where A = 27r/k is the electron wavelength. The result for the A 3> -R^ case is vahd in the 
Born approximation for scattering by well-separated impurities; ^'^ in the opposite limit the 
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electron is moving semiclassically in a smooth potential, where r oc depends on the 

electron energy. 

To describe the magnetic field, we choose the Landau gauge A = (0, Bx, 0) where the 
eigenstates \nKs) of the unperturbed Hamiltonian p|/2m + /x^gf (cr • B) /2 are represented 
by the spinors: 

where fig is the spinor corresponding to one of the states |s) = |T)Ji)) ^ is the Landau level 
number, the magnetic length Ib = ^hc/\e\B, xk = x — Xk, Xk = —kyl% is the center of 
the oscillator wavefunction, ky = ~2ttK / Ly, = 0, 1, . . . , -f^max, Ly is the axis size of the 
system and Hn is the nth Hcrmitc polynomial. The corresponding unperturbed spectrum is 
E{n, s) = hwc{n + 1/2) ± giisB /2. Due to the selection rules for the matrix elements of py 
and spin components (Tx and cr^, the SO coupling only connects states of opposite spins from 
nearest Landau levels. This results in a small shift in the energies of the order oim{a/h) . 
The matrix elements of the disorder Hamiltonian diagonal over the spin index are given by: 

H,^^{n'K's';nKs) = J d^p^,,j,,Ap)U{p)<PnKs{p), (5) 

and couple states in all Landau levels for which \Xk — Xk'\ < Ib, and, correspondingly 
\ky — k'y\ < thus leading to electron localization both in the x— and y— directions. As 
a result the density of eigenstates with energies Ej has the form of broadened Landau levels, 
typically {at B = 5 T) with the width F on the order of or less than 1 meV, at least an order 
of magnitude smaller than hcuc- The corresponding eigenfunctions are expressed as linear 
combinations: 



H„ 
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1^.) = E «nL 1^^^) , (6) 
n,K,s 

with complex coefficients a^^Ks- Only a few of the amplitudes a^^^ are nonnegligible, since 
U (p) couples only the states with close momenta ky. Due to the combined effect of disorder 
and SO coupling, the \nKi |) and 1^X2 i) states within one Landau level become coupled 
as shown in Fig. 1(a), thus introducing randomness in the spin precession at the frequency 
scale T/h, and, in turn, spin relaxation. 

We consider an electron-hole plasma injected by a light pulse in an undoped quantum 
well (Fig. 1(b)). We assume that the plasma density is small, and neglect all many-body 
effects. Alternatively, for a doped quantum well, we assume the carriers are injected into 
unoccupied Landau states above the Fermi level. We concentrate on the spin relaxation of 
electrons, since the hole spins relax much faster. 

We assume that the spectral width of the exciting light Auj satisfies the conditions 
F <^ hAuj <^ fkUc as shown in Fig. 1(b). Therefore, the states \ip"^{t = 0)) (index 1 < i < N^n 
numerates the injected electrons, with their total number being N[n) in which the elec- 
trons are injected can be written as superposition of the \nKis) states from essentially one 
Landau level and spin projection determined by the light polarization.^^ Within 
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this approximation a unitary transformation connects the full sets of \iljf^{t = 0)) and 
IniiiA'jSin) = \^kX^ ~ 0)) states. Thus, in place of averaging the time-dependent spin com- 
ponents over \ipl'^{t)), we can average over the as follows: 

('^cW) = TrE(V'rwkclCW> (7) 

i=l 

' in ■ 1 

1=1 

This relaxation is coherent in the sense that the total energy (H) of the injected electrons 
is conserved on the timescale considered here and is not transformed into lattice phonons or 
low-energy electron excitations. 

To evaluate numerically, we have to find the spinor representation of the 1$"^. (t)); 

that is, we must determine the eigenstatcs of the Hamiltonian H in Eq.(l). To do 
this, we have chosen a basis of = 256 states per Landau level for each spin pro- 
jection. The total Hilbert space included — 6 Landau levels, which we find suffi- 
cient for our choice of the parameters given below, for which HlOc ^ T- The coefficients 
^■liKs considered as the elements of the eigenvectors are arranged in the following order: 
({oo/ct) ' {^oKi} {(^Nl-iki} , {<^Nl-iki}) , where in every subset K is running from to 
K 



III. NUMERICAL RESULTS: THE ROLE OF THE (/-FACTOR AND THE DIS- 
ORDER LANDSCAPE 

For sample calculation we consider two types of typical structures. The first is a symmet- 
ric quantum well with electrons located in a GaAs layer, with the Dresselhaus SO coupling 
ckd = 0.35 X 10~^ eVcm, m — 0.067mo (mo is the free electron mass), and g — —0.45. The 
other is an asymmetric structure with the electrons located in a Ino.5Gao.5As layer, with^^ 
a Rashba SO coupling cir = 0.35 x 10~^ eVcm, m = 0.05mo, and (7 = 4. We assume that 
both of them have a width of 10 nm. For these structures the ratio of effective masses 
m(GaAs)/m(lno.5Gao.5As) = 1.35, while the electron factors arc different approximately 
by a factor of |(7(lno.5Gao.5As)/5'(GaAs)| = 8.9. It is the difference in (^-factors that will lead 
to large quantitative difference between the spin relaxation in these two structures. 

As a realization of the disorder^^ we consider a short-range potential with Vg — 3.5 meV 
for GaAs and Vg — 5.0 meV for Ino.5Gao.5As structure, Rg = 3 nm, and iVimp = 10^^ cm~^. 
These parameters can describe the random potential arising due to the single unit cell layer 
thickness variations in quantum wells of the width of 10 nm. Since Ib ^ Rg, the width of 
the Landau level F = U^IlOcIt^t does not depend on the level number. We also consider a 
long-range potential with Rg — 20 nm and the same amplitudes Vg as for the short-range 
one. Both these random potentials lead to mobilities of 5 x 10^ cm^/Vs at concentrations 
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of electrons A^ei ~ 5 x 10^^ cm~^ in both structures. Note that in the case of a long-range 
potential Rg ^ Ib, we have T = 2^/JJf^. This spatial scale can be achieved by quantum 
well fabrication with growth interruption^^, or by remote doping on the quantum well sides. 

The results of numerical calculations for the GaAs quantum well are presented in Figs. 2 
and 3 for different initial spin components {(TziGi)) — —1 and ((Ta;(0)) = 1, respectively, two 
different types of disorder, and magnetic fields of 4 T and 8 T. In these GaAs-based structures 
we find the relaxation times on the order of 100 ps, of the same order of magnitude for out-of- 
plane {az{t)) and in-plane {(7x{t)) spin components. They relax on the same time scale since 
the same mechanism, that is the random spin precession, leads to the relaxation in both spin 
polarizations. With the increase of the Landau level number, the electron motion becomes 
less sensitive to the disorder, and the spin relaxation time increases, as demonstrated by the 
results for n = 2 in Fig. 2 and Fig. 3. This corresponds well to the experimental results of 
Sih et al?^ Spin relaxation due to acoustic phonon emission^^ occurs at a much longer time 
scale, and is not considered here. An interesting effect is the spin precession in the {(Tz{t)) 
relaxation, clearly seen in our calculations. The spin-orbit interaction couples Zeeman-split 
Landau levels which are broadened due to disorder, and so the ;z-component of the spin is 
not a constant of motion. For this reason, the initially prepared ci^ = ±1 states precess with 
the frequencies in the range determined by the Zeeman splitting and the Landau level width, 
distributed over interval of the width V/h centered at qii-qB /h. The precession amplitude is 
smeared with time due to spin relaxation, as seen in our results. 

Fig. 4 presents the results for Ino.5Gao.5As structure, where the spin relaxation is much 
slower. Due to the increased Zeeman splitting the frequency of the spin oscillations here 
is considerably larger than that in the GaAs quantum well, and the oscillations become 
more well-defined. The role of the Zeeman splitting \g\^BB in spin relaxation, clearly 
seen when one compares the results for the GaAs and lno.5Gao.5As structures, is crucial 
to the mechanism of the spin relaxation, and can be understood as follows. The random 
potential accompanied by SO coupling leads to the spin-flip transitions. However, only 
spatially close states (with a large overlap) with opposite spins can contribute effectively to 
the spin relaxation. On the other hand, the spin-fiip process should conserve energy, and, 
therefore, a lateral distance 4 on which the orbital has to be displaced to find its spin-fiip 
partner state depends on the Zeeman splitting. To understand the effect we evaluate the 
fluctuation of the expectation value of energy for a state described by wave function -0 (p) : 



((AC/p)^> 



(8) 



C/(pi) (V'' (Pi) - V'' (Pi + p)) ciVi 

= 2{U^) J F,{p, - p2)^' (pi) [V'' (P2) - V'' (P2 + P)] d^Pid'p2. 
This fiuctuation can be calculated in two limiting cases as: 

{{AU,f)^2{U'){ ^2 , _ (9) 



I'B '- -' 
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where we assumed for an estimate that for the ground orbital state ip"^ (p) = 
exp (— /27r/^. If 4 determined by energy conservation is less than Ib, that is 
\g\fiBB < {U'^)^^'^ mhi {I B / Rg, Rg/ Is) , the relaxation occurs effectively. In the opposite 
case, the decay time increases due to a small spatial overlap of the initial and final states. In 
other words, if \g\iXBB > T, the spin relaxation is suppressed, since the energy conservation 
cannot be fulfilled in the spin-flip process. In this case the 2;-component of spin relaxes from 
the initial value by approximately r/\g\iiBB, and then the system has to pass by emitting 
acoustic phonons through a phonon bottleneck^^ for the spin to relax it further. 

This "spin-fiip distance" argument explains the difference between the spin relaxation for 
a "long" and "short"- range potentials shown in Fig. 4. In both cases the amplitude of the 
potential fluctuations is the same, but in the case of a long-range potential the electron must 
be displaced a longer distance to find its spin-flip partner state. Thus, the relaxation rate 
decreases with the increase of the correlation length Rg, in agreement with results presented 
in Fig.4. 

IV. CONCLUSIONS AND POSSIBLE APPLICATIONS 

To conclude, we investigated the spin relaxation in quantizing magnetic fields in a disor- 
dered 2D electron gas, and found that the result depends on the Landau level number, and 
that the process is accompanied by spin oscillations. The mechanism discussed here is dif- 
ferent from the usual Dyakonov-Perel' mechanism, and is closer in nature to the Elliot- Yafet 
mechanism, in which the spin relaxation rate increases with the disorder.^ The relaxation 
rate depends on the details of the potential and electron gi— factor, and cannot be understood 
solely in terms of the electron mobility. 

The results obtained can be used to characterize the disorder in undoped quantum wells. 
Such disorder is typically probed experimentally by studying the inhomogeneous broadening 
of the spectra of cxcitons, a technique restricted either to the spatial scale given by the 
exciton Bohr radius or by the exciton localization length. Neither of these lengths can be 
changed externally in a well-controlled way. The advantage of studying the spin relaxation 
in a magnetic field is that it allows probing different spatial scales of disorder by controUably 
changing the length Ib trough varying the magnetic field. 
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FIG. 1: (a) Schematic plot of the broadened Landau levels and interlevel transitions for or = 
0, ao 7^ 0. Arrows label electron spins, and the index n corresponds to the Landau level number, 
(b) Schematic plot of the optical transitions between the Landau levels. Only one spin projection 
is presented. 
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FIG. 2: (Color online) (a) (cr^(t)) for a short-range {Rg = 3 nm) potential, and (b) {(Jz{t)) for a 
long-range {Rg = 20 nm) potential in the GaAs quantum well. Solid lines correspond to n = 0; 
the dashed line in Fig. 2(a) corresponds to n = 2, B = 4T. The strength of magnetic field is shown 
near the lines. 
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FIG. 3: (Color online) {(7x{t)) for a short-range potential {Rg = 3 nm) in the GaAs quantum well, 
with B = 4 T. Landau level numbers are shown near the lines. 
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FIG. 4: (Color online) {(Jz{t)) for long-range and short-range potentials in the Ino.sGao.sAs quan- 
tum well. 
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